Effect of Strong Electron Correlation on the Efficiency of Photosynthetic Light 

Harvesting 
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Research into the efficiency of photosynthetic light harvesting has focused on two factors: (1) 
entanglement of chromophores, and (2) environmental noise. In this Letter efficiency is shown 
to be equally affected by a third factor: (3) strong electron correlation. While chromophores are 
conjugated 7r-bonding molecules with strongly correlated electrons, previous models have not treated 
this correlation explicitly. Here we generalize the single-electron models to a multi-electron model 
and show that correlation enhances the energy-transfer efficiency in the model by more than 100%. 
Implications include insights into the interplay of correlation and entanglement as well as new design 
principles for energy-efficient materials. 

PACS numbers: 31.10.+Z 



Nature harvests solar energy with a remarkably high 
quantum efficiency, the percentage of charge carriers cre- 
ated by photons. Recent spectroscopic experiments [D- 
3] and theoretical models |I-13|, provide strong evi- 
dence that efficient light harvesting in nature occurs by 
a quantum mechanism involving sustained electronic co- 
herence 14 1 and entanglement [l6| between chro- 
mophores. While the chromophores are chlorophyll 
molecules containing large networks of conjugated car- 
bon bonds that surround a charged magnesium ion, 
they have largely been represented in theoretical stud- 
ies [4j4l3j by one-electron models that neglect the effects 
of electron correlation and entanglement within chro- 
mophores. Two advanced methods in electronic struc- 




ture, density-matrix renormalization group 17| and two- 
electron reduced-density-matrix theory IH l3> have re- 
cently shown that networks of conjugated bonds as in 
acene chains [l7|, H3| ; acene sheets [20| , and chlorophyll 
are associated with polyradical character that cannot 
be adequately described without a strongly correlated 
many-electron quantum model. 

In this Letter we examine the efficiency of light har- 
vesting where we represent each chromophore by a corre- 
lated TV-electron model to treat strong electron correla- 
tion. Figure 1 illustrates the replacement of one-electron 
models for each of the 7 chromophores in the Fenna- 
Matthcws-Olson (FMO) complex of green-sulfur bacte- 
ria by TV-electron models of increasing complexity. Here 
we em ploy the TV-electron Lipkin-Meshkov-Glick (LMG) 
model [21| used extensively in electronic structure 2^- 29j 
and quantum information [l6| . We find that strong elec- 
tron correlation dramatically enhances the efficiency of 
the energy transfer to the reaction center by more than 
100%. Strong electron correlation, the results show, is 
likely employed by nature to enhance its energy-transfer 
efficiency as much as other factors such as (i) environ- 
mental noise 1 5|— KL 0I| and (ii) entanglement between chro- 
mophores llMl3| , which have been extensively studied in 
the recent literature. The correlation of electrons within 



FIG. 1. Single electrons or correlated chromophores. 

Each of the seven chromophores in the FMO complex is gen- 
erally treated as a single electron in a two-state model (left), 
and yet the chromophores are constructed from chlorophyll 
molecules with many strongly correlated electrons (right). 
Here we treat each of the chromophores by a correlated N- 
electron model by Lipkin, Meshkov, and Glick. Illustration by 
K. Naftchi-Ardebili, The University of Chicago, 2011. Used 
with permission. 



a molecular subunit like a chromophore cannot be sep- 
arated from their entanglement between molecular sub- 
units. Furthermore, the results suggest a general design 
principle for man-made materials in which electron corre- 
lations and entanglements both within and between sub- 
units are simultaneously tuned for achieving enhanced 
quantum efficiencies. In combination with other recent 
advances, including the study of functional subsystems 
of the FMO complex this interplay of electron en- 
tanglements on different length scales may enable us to 
develop materials with quantum efficiencies approaching 
those found in natural processes from photosynthesis to 
bioluminescence. 

The Fenna-Matthews-Olson (FMO) complex of green- 
sulfur bacteria contains three identical subunits, each 
with a network of seven chromophores. Theoretical mod- 
els of an FMO complex's subunit typically represent each 
of the seven chromophore by a one-electron model in 
which the electron has access to two energy levels sep- 
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arated by the excitation energy of the chromophore. In- 
teractions U between pairs of chromophores are modeled 
by the exchange of single-electron excitations (or exci- 
tons) between them: 



H = H + U 



where 



(i) 



(2) 



U = ^2u S}t al +1 a s ^ 1 al_ 1 a ty+ i. (3) 

The quantum numbers s and t denote the seven sites of 
the chromophores while the quantum number to, equal 
to +1 or -1, indicates one of the two energy levels within 
each chromophore. The second-quantized operator a\ m 
(a Sjm ) creates (annihilates) an electron on chromophore 
s in energy level to. The 7 parameters e s are the excita- 
tion energies of the chromophore, and the 21 parameters 
U s> t are the coupling energies between all pairs of chro- 
mophores. Typical values for the excitation and coupling 
energies are given in the 7x7 Hamiltonian of Ref. 30] . 

Here we generalize this representation by using the N- 
electron LMG model [2l[ . Each chromophore is modeled 
as N electrons in two energy levels that are each A^-fold 
degenerate; a pair wise interaction scatters two electrons 
from the lower level to the upper level or from the upper 
level to the lower level 24J. The Hamiltonian of the 7 
interacting LMG chromophores can be expressed as fol- 
lows: 



H = Hn + XU + V 



where 



H o = ~ me s a t,m, P a s 



(4) 



(5) 



U = 



X/ ^M a l, + l,p a S,-l:P a t,-l, P a t, + l,P' (6) 



\,m,p a \.m,q a s.-m,qO's, — m,p-, 



(7) 



s,m,p,q 



where V accounts for the interactions within each chro- 
mophore and A G [0, 1] is a screening parameter to be 
defined below. The new quantum number p (pr q) de- 
notes the N degenerate states within each energy level, 
which are necessary to accommodate the N electrons. 
When the interaction strength V equals zero, the N elec- 
trons on each site are non-interacting, and the model be- 
haves the same as the one-electron-per-site models; when 
V is unequal to zero, we have a generalized model for 
light harvesting with a tunable electron correlation on 
the chromophores. 

Environmental effects of dephasing and dissipation as 
well as the transfer of energy to the reaction center (sink) 
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FIG. 2. Populations of chromophores 1-3 and sink with 
(b) and without (a) electron correlation per site. The 

exciton populations in chromophores 1, 2, and 3 as well as 
the sink population are shown as functions of time in fem- 
toseconds (fs) for N = 4 with (a) V = 0.0 and (b) V = 0.8. 
Correlating the four electrons on each chromophore signifi- 
cantly accelerates the increase in the sink population with 
time. 



can be incorporated by introducing a Lindblad operator 
L into the quantum Liouville equation 



dt 



L(D) 



(8) 



where D is the many-electron density matrix. The Lind- 
blad operator can be divided into three operators that 
account for dephasing, dissipation, and the sink 



L{D) = L deph (D) + L diss (D) + L sink {D) (9) 



where 



L deph (D) = aJ2^(k\D\k)\k}(k\-{\k)(k\,D}, (10) 

k 

L diss (D) = pJ2 2 ( k \ D \ k )\9)(9\ - {\k)(k\,D}, (11) 

k 

L sink (D) =27(3|D|3)|«)(s|-7{|3)(3|,D}. (12) 

The state \g) denotes the ground eigenstate of the Hamil- 
tonian in Eq. ((4]), the states \k) represent the excited 
eigenstates computed from this Hamiltonian where the 
interaction U between chromophores is set to zero, the 
state \s) denotes the reaction center (sink), and |3) indi- 
cates the first excited state of the third chromophore mul- 
tiplied by the ground states of the other chromophores. 
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Non-Markovian effects from the environment can also be 
added to the quantum Liouville equation (for example, 
refer to Ref. [HI); however, they will not qualitatively 
change the effect of strong electron correlation within 
chromophores on energy-transfer efficiency. 

The values for the site excitation energies are defined 
by the 7x7 Hamiltonian of Ref. [30] ; when V > 0, the 
parameters {e s } are adjusted to ensure that the excita- 
tions of the correlated LMG models agree with the exci- 
tation energies of the seven chromophores. Selection of 
the coupling energies U s>t in Eq. 2]) also requires adjust- 
ment. Because the coupling energies given in Ref. [13] 
are "dressed" dipole-dipole interactions that account im- 
plicitly for both the electron correlation within the chro- 
mophores and the protein environment surrounding the 
chromophores, they require adjustment for the LMG 
chromophore model that contains an explicit electron- 
electron interaction. Previous 2-RDM-based calculations 
of acene chains and sheets reveal the presence of strong 
electron correlation in conjugated systems of sizes similar 
to the bacteriochlorophyll molecules in the FMO com- 
plex. Based on acene-chain data [2(J, we estimate the 
interaction strength V with N = 4 to be 0.8, which gives 
an 80% probability of finding an electron in the highest 
occupied orbital (occupied in a mean-field treatment) and 
a 20% probability of finding an electron in the lowest un- 
occupied orbital. This estimate is conservative because: 
(i) the acene chains of a similar length typically reveal 
a nearly biradical filling (« 50% in the highest occupied 
orbital), and (ii) the presence of the Mg ion with its d 
orbitals is expected to enhance the degree of correlation. 
To prevent overcounting of the electron correction's ef- 
fect on the coupling, we screen the coupling energies C/ S)t 
of Ref. 30] by selecting A in Eq. [4j to be less than unity. 
Specifically, we set A = 0.629 to match experimental data 
with the LMG model when N = 4 and V = 0.8. The rate 
parameters a, /3, and 7 in the Lindblad operators in Eq.[9] 
are chosen in atomic units to be 1.52xl0~ 4 , 7.26xl0 -5 , 
and 1.2fxf 0~ 8 , respectively. These definitions are similar 
to those employed in previous work when V = [30] . 

The exciton populations in chromophores 1 , 2, and 3 
as well as the sink population are shown as functions of 
time in femtoseconds (fs) in Figs. 2a and 2b for N = 4 
and A = 0.629 with V = 0.0 and V = 0.8, respectively. 
Population dynamics of the excitation are generated by 
evolving the Liouville equation in Eq. (|SJ) from an initial 
density matrix with chromophore 1 in its first excited 
state and the other chromophores in their ground states. 
Correlating the 4 electrons on each chromophore signifi- 
cantly accelerates the increase in the sink population with 
time. By 1 ps the sink population for V — 0.0 is 0.114 
while the population for V — 0.8 is 0.287. Correlating 
the excitons on each chromophore also has the effect of 
shortening the periods of oscillation in chromophores 1 
and 2 and accelerating the population decay in these 
chromophores, which is consistent with the change in the 



sink population. 

The sink population as a function of time (fs) is shown 
in Figs. 3 for a range of V with N = 4 and A = 0.629. 
Importantly, as V increases, we observe a dramatic ac- 
celeration of the increase of the sink population. For 
V increasing by the sequence 0.0, 0.4, 0.8, and 1.2, the 
sink population at 2 ps increases by the sequence 0.221, 
0.367, 0.498, and 0.547. For N = 4 correlating the elec- 
trons within the chromophores significantly increases the 
efficiency of energy transfer to the reaction center (sink) 
by as much as 148%. While we choose N = 4, the num- 
ber TV of electrons per chromophore can model electron 
correlation for any N > 2. The precise value of N > 2 
is unimportant because the effect of changing N can be 
related to a rescaling of the interaction V. 
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FIG. 3. Correlation-enhanced transfer to the reaction 
center. The reaction center (sink) population as a function of 
time (fs) is shown in for a range of V with N = 4. Correlating 
the electrons within the chromophores significantly increases 
the efficiency of energy transfer to the reaction center (sink). 

Figure 4 examines the entanglement of excitons be- 
tween the LMG-model chromophores for N = 4 and A = 
0.629 with V = 0.0 and V = 0.8. We employ a measure 
of global entanglement in which the squared Euclidean 
distance between the quantum density matrix and its 
nearest classical density matrix is computed [2(| [3l| - l33l | : 



a(D) = \\D -C\\ 2 =J2 



(13) 



In some cases like the entanglement of the chromophores, 
the squared Euclidean distance can be viewed as the sum 
of the squares of the concurrences 15| , a measure of local 
entanglement. Within the mathematical framework of 
Bergmann distances, the squared Euclidean distance can 
also be related to quantum relative entropy 13j,ll6|, which 
is often applied as a global entanglement measure. The 
squared Euclidean distance er(Z?) is nonzero if and only 
if the excitons on different chromophores are entangled. 
The correlation of electrons increases the degree of the 
entanglement between chromophores at early times and 
the frequency of its oscillation. The greater entanglement 



4 




200 
Time (fs) 



300 



400 



FIG. 4. Entanglement of excitons with and without 
electron correlation. A measure of global entanglement is 
shown as a function of time (fs) for N = 4 with V = 0.0 and 
V = 0.8. The correlation of the excitons increases the degree 
of the entanglement between chromophores at early times and 
the frequency of its oscillation. 



at early times reflects the opening of additional channels 
between chromophores for quantum energy transfer. 

The chromophores interact through intermolecular 
forces, both dipole-dipole and London dispersion forces. 
The present results imply that nature enhances these in- 
termolecular forces through strong electron correlation in 
the 7r-bonded networks of the chromophores to achieve 
the observed energy-transfer efficiency. The two parame- 
ters of the LMG chromophore model provide the simplest 
approach to studying the effect of strong electron correla- 
tion V on the effective coupling between chromophores. 
The one-electron or dipole models with their coupling 
energies U s ,t can mimic the efficiency from such correla- 
tion within chromophores through an empirical inflation 
of the one-electron U coupling, but they do not provide a 
mechanism for either isolating or estimating the magni- 
tude of the enhanced coupling due to strong electron cor- 
relation. Orbital occupations from recent 2-RDM calcu- 
lations of correlation in polyaromatic hydrocarbons (20j 
suggest V — 0.8 to be a conservative estimate of the cor- 
relation within the LMG models of the 7 chromophores. 
Using this estimate with coupling energies screened to 
match experimental and computational data, we observe 
a greater than 100% enhancement from the strong elec- 
tron correlation. 

Correlation-enhanced energy transfer can be compared 
with noise-assisted transfer. Theoretical models 0t3,I3 
[ili ] have shown that noise from the environment (de- 
phasing) can assist energy transfer in the FMO complex 
by interfering with the coherence (resonance) between 
chromophores with similar energies, which facilitates the 
downhill flow of energy to the lowest-energy, third chro- 
mophore, connected to the reaction center. Electron 
correlation on each chromophore, we have shown, en- 
hances transfer by opening additional coherent chan- 
nels between chromophores, which also accelerates en- 



ergy transfer to the third chromophore. Photosynthesis 
can draw from both of these sources, strong electron cor- 
relation within chromophores and environmental noise, 
to increase the rates of energy transfer. Some of the ex- 
perimental energy-transfer efficiency attributed to noise 
in one-electron chromophores models may in fact be due 
to strong electron correlation. 
Briggs and Eisfeld 



12| recently examined whether 



the energy-transfer efficiency from quantum entangle- 
ment might be matched by a purely classical process. 
They conclude that if chromophores are approximated 
as dipoles, then quantum and classical treatments can 
achieve similar efficiency. While their result might also 
be extendable to other dipole or one-electron approxi- 
mations of the chromophores, many-electron models of 
the chromophores cannot be represented within classical 
physics. Neither the electron correlation, present in the 
LMG model of the chromophores when V > 0, nor the 
associated enhancement of energy-transfer efficiency can 
be mapped onto an analogous classical process. 

Many theoretical models (BpS 12, 13 1 have been de- 
signed to explore the energy transfer in recent light- 
harvesting experiments, but most of them treat each 
chromophore by a single electron with two possible en- 
ergy states. In reality, however, the chromophores are 
assembled from chlorophyll molecules that contain an ex- 
tensive network of conjugated carbon-carbon bonds sur- 
rounding magnesium ions, from which significant strong 
electron correlation, includi ng p olyradical character, has 
been shown to emerge 17J, |20( . In this Letter we have 



examined the effect of strong electron correlation and en- 
tanglement within chromophores through an extension of 
single-electron models of the chromophores to iV-electron 
models, based on the LMG model [13, El, |2fj| . We 
find that increasing the degree of electron correlation of 
each LMG-model chromophore significantly enhances the 
efficiency with which energy is transferred to the reaction 
center (sink). This result, showing that nature likely uses 
strong electron correlation to achieve its energy-transfer 
efficiency, also has implications for the design of more 
energy- and information-efficient materials. 
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